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Abstract 

Contributions of D waves to physical observables for neutral pion photoproduction from the 
proton in the near-threshold region are studied and means to isolate them are proposed. Various 
approaches to describe the multipoles are employed — a phenomenological one, a unitary one, 
and heavy baryon chiral perturbation theory. The results of these approaches are compared and 
found to yield essentially the same answers. D waves are seen to enter together with S waves 
in a way that any means which attempt to obtain the Eq + multipole accurately must rely on 
knowledge of D waves and that consequently the latter cannot be dismissed in analyses of low- 
energy pion photoproduction. It is shown that D waves have a significant impact on double- 
polarization observables that can be measured. This importance of D waves is due to the soft 
nature of the S wave and is a direct consequence of chiral symmetry and the Nambu-Goldstone 
nature of the pion. F-w&ve contributions are shown to be negligible in the near-threshold region. 
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I. INTRODUCTION 



Due to the spontaneous breaking of chiral symmetry in Quantum Chromodynamics 

n 

(QCD) the 7r meson appears as a pseudoscalar Nambu-Goldstone boson As a dynamical 
consequence, the S"-wave amplitude for the — > tt°N reaction is small in the threshold 
region, because it vanishes in the chiral limit, i.e. when the light quark masses are set equal 
to zero . An additional consequence is that the P-wave amplitude is large and leads 

to the A resonance at intermediate energies [5| . Accordingly, the photoproduction of neutral 
pions differs from the general pattern for hadronic reactions where the S wave dominates 
close to threshold and then, as the energy increases, the higher angular momentum waves 
(P, D, . . . ) start to become important. By contrast, for the 7 A" — > 7r°N reaction the S- and 
P-wave contributions are comparable even very close to threshold {(J. Hence, the accurate 
extraction of S and P waves from pion photoproduction data becomes an important issue 
in the study of the breaking of chiral symmetry and QCD. 

The partial waves (electromagnetic multipoles) are not experimental observables, but 
rather are quantities extracted from data using some kind of approach and/or theoretical 
input — Heavy Baryon Chiral Perturbation Theory (HBCHPT), for instance. 

'. n the first resonance region (A region) the influence of higher partial waves was addressed 
in ?J employing the energy-dependent solution of SAID 8]; however when it comes to 
precision physics in the near-threshold region, the standard approach consists in analyzing 
data assuming that the partial waves of interest (S and P waves) are sufficient to describe 

the experimental data so, higher partial waves may be neglected q, y, y, y, y, q. ^ 

the literature there are several phenomenological models of neutral pion photoproduction 
based on the standard Born terms [ll| that are applied at energies spanning from threshold 
to the first resonance region. These models compute the electromagnetic multipoles and 
include higher partial waves in their analysis, although they are not adequate to address 
precision physics in the near-threshold region due to model dependencies (e.g. form factors 
and resonances treatment) and large uncertainties in the multipoles. As an example, the 
most accurate model calculation in the near-threshold energy region is the Dubna-Mainz- 
Taipei (DMT) dynamical model 12| whose P waves do not provide a good description of 
the experimental data for the polarization asymmetry. Two arguments have motivated this 
assumption in previous analyses: (i) The angular dependence of the experimental differential 
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cross section can be described accurately using Legendre polynomials up only to order two. 
Indeed, S and P waves constitute the minimal set of partial waves needed to reach that 
angular dependence; (ii) In the near-threshold energy region higher partial waves are weak 
and the early dominance of the M 1+ multipole renders them negligible. 

As will be discussed in Sec. |TTJ the first statement is misleading because any coefficient 
that accompanies a Legendre polynomial in the expansion depends on all partial waves, im- 
plying that in some circumstances higher partial waves can make an important contribution 
to the coefficients that accompany the lower-order Legendre polynomials, thereby posing 
uncertainties in the multipole extraction. This assumption cannot be taken for granted and 
has to be tested. The second argument is actually in favor of the possible importance of 
higher partial waves and not against it, as is usually stated. The conclusion disregards that 
the S wave is also weak and that important features of the angular dependence of the ob- 
servables are dominated by the interference of different partial waves. In this situation the 
dominance of a certain contribution (such the M\ + in this case) can lead to an important 
enhancement of smaller partial waves through interference, making them relevant. An ex- 
cellent example of this situation is the well-known S-w&ve/ D-w&ve interference in the beam 
asymmetry in 7] photoproduction, where a very weak D wave turns out to be necessary to 
explain the experimental data due to its interference with the dominant S wave. In this 
article it is shown that something similar, although more subtle, happens in the case of pion 
photoproduction, where D waves affect the extraction of the S wave. 

The purpose of this article is twofold: 1) to complement and extend the analysis of what 
was reported in Ref. [k^] where the impact of .D-wave contributions in the near-threshold 
region was explored using as starting point HBCHPT and 2) to study the impact of the 
.D-wave contribution on the observables and assess the measurability of this impact. 

As will be seen in next sections the understanding of D waves and their impact on the 
observables is necessary if one hopes to extract the S wave accurately and to arrive at 
conclusions regarding isospin breaking, unitarity, energy dependence of the multipole, the 
magnitude of the unitary cusp, or the extraction of the Low Energy Constants (LECs) in 
HBCHPT. 
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II. GENERAL STRUCTURE OF THE OBSERVABLES IN TERMS OF THE 
ELECTROMAGNETIC MULTIPOLES 



It has been customary to study pion photoproduction from the nucleon in the near thresh- 
old energy region (up to approximately 170 MeV) using only the contributions of S and P 
waves and there is abundant literature that discusses the structure of the observables in 
terms of only these partial waves —namely, using E 0+} Ei + , Mi_, and 

M 1+ electromagnetic multipoles. In Ref. 13| it was found that D waves play an unex- 
pected role in this energy region and that D waves are required to extract the Eq + multipole 
accurately. 

This article expands the work reported in Ref. [3| and the first step consists in analyzing 
the structure of the observables in terms of the reported electromagnetic multipoles when the 
D waves are added. All photoproduction observables can be written in terms of responses of 
two kinds, time reversal even (TRE) and time reversal odd (TRE) . Both types of responses 
are real quantities obtained by taking either the real or imaginary part of bilinear products 
of the complex multipoles, respectively. If the real part is chosen, the observable is TRE 
(cry, E, E, F, C x , C z , T x , T z , L x , and L z ) and if the imaginary part is chosen the observable 
is TRO (T, P, G, H, O x , O z ). See Refs. 15|, [16] for the definition of the observables. From 



this full set of observables, in this article the focus is placed on those that do not require 
one to measure the recoil polarization, in other words, those that involve polarizing only the 
beam, the target, both or neither, namely the ones that could be accessed experimentally 
in the near-threshold region. That leaves four TRE — cry, S, E, and F — and four TRO — T, 
P, G, and H — observables, which will be discussed in more detail in the following. In Sec. 
IIVI a selection of several observables at different energies and angles is presented, restricting 
the discussion (and the tables in Appendix |A|) to the most promising ones involved in the 
search for D waves and their interplay with S and P waves, namely, the differential cross 
section <tt, the photon beam asymmetry E, and the double-polarization asymmetries E and 
F (circular photon polarization and polarization of the target, respectively, along the beam 
axis and the orthogonal axis within the scattering plane). For the nomenclature used see 
Table [U the structure of these observables is provided in Appendix El All these observables 
have become accessible experimentally in recent years thanks to the existence of high-duty- 
cycle photon facilities such as MAMI (Mainz) and HI7S (Duke) and to the development of 
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TABLE I: Relationships among several notations present in the literature for the asymmetries 
and responses. In this article we use the sign conventions of Ref. lf| for the asymmetries. The 



first line contains different notations for the differential cross section and the associated response. 
The second line stands for the photon beam asymmetry, measured as the difference between the 
differential cross section for perpendicular and parallel polarized photons. The third line contains 
several conventions for the beam-target asymmetry obtained by measuring the difference of the 
cross sections for the two circular polarizations of the photon with a target polarized along the beam 
direction. The fourth line contains several conventions for the beam-target asymmetry obtained 
by measuring the difference of the cross sections for the two circular polarizations of the photon, 
but now with target polarized in the sideways direction. In the second column, 7 stands for linear 
polarization of the photon while j c stands for circular polarization of the photon. The z stands 
for longitudinal polarization of the target (along the photon direction) and the x for sideways 
polarization of the target (transverse and in the scattering plane). In the third column, in W T (0) 
and W TT (0) stands for unpolarized target. The s in W T> (s) stands for target polarized in the 
sideways direction (in any other place in this article s stands for the Mandelstam variable) and the 
/ in W T ' (I) for target polarized in the beam direction. 

This article Ref. [10] Ref. [14] Ref. [15] 

or = qirWr/kj ; Eq. {SI) a T = q^Rf /ky q 7r W T (0)/k 1 da/dfl 

£ = -Ws/Wt ; Eq. flUD Atf) = -nSgr/RSp -W TT (0)/W T (0) £ 

E = We/W t ■ Eq. (H3} A{ lc ,z) = -R° T Z T ,/R™ -W T ' (I) /W T (0) E 

F = W F /W T ; Eq. flM} -A( lc ,x) = -R^ X T ,/R^ W r (s) /W T '(0) —F 



polarized targets. 

Any pion photoproduction response 71 (s, 6) depending on Mandelstam variable s and 
photon-pion angle 9 can be expanded in terms of Legendre polynomials Pj (8) times sin 9 to 
a specific integer power n: 

n(s,9) = sm n 9[TZ (s) +TZ 1 (s)V 1 (9) + . . .} , (1) 

where coefficients IZj (s) can be denned on terms of the electromagnetic multipoles up to a 
certain partial wave and the number of Legendre polynomials also depends on up to which 
partial wave the observable is expanded. One has n = for Wt and We, n — 1 for Wp and 
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n = 2 for Ws- 

For example, the differential cross section can be written including up to D waves as: 
a T (s, 9) = ^ [T (s) + 7\ (s) V 1 (9) + T 2 (s) V 2 (9) 

+T 3 (s)V 3 (9)+T A (s)V 4 (9)} , 

where q n and fc 7 are the pion and photon momenta in the center of mass, respectively, and 
the Tj depend on the photon energy. 



The full calculation of the responses up to D waves can be found in Ref. [17|, although 
here we have preferred to proceed using the developments in Ref. 14j and to present the 
results using tables instead of lengthly equations (see Appendix [X]) . The reason for using 
tables will become apparent later when it will be seen that they make it easier to extract 
the interplay of the interferences between partial waves. 

The first place where effects of D waves in low-energy pion photoproduction were found 
[kJ ] was the part of the differential cross section that is associated with the Legendre poly- 
nomial V\ (9), namely T\ in Eq. (T2]). This is used here as a working example of how to 
employ the tables. From Eq. (1A9I) this piece of the response can be written: 



Ttis) = J^Re{M*(s) T? M^s)} , (3) 

where M.j (s) = E 0+ , E 1+ , E 2+ , E 2 -, Mi_, M 2+ , M 2 _, and the coefficients T[ 3 can be 

read from Table IIVI in Appendix |A], obtaining 

T 1 = 2Re[P*E 0+ }+6T 1 , (4) 

where P\ = 3E± + + M\ + — Mi_, and 5T\ stands for the D-wave/P-wave interference contri- 
bution 

27 



8T X = 2 Re 



.-.U,' . .U,. + (M* + - M*_) /•;, 



f72 3 9 9 

+ E* 1+ ^-E 2+ - -P 2 _ + -M 2+ - -M 2 _ J ( 5 ) 

+ 0M* + + 3M;„^) M 2 

Proceeding in a similar way, the other observables and multipolar expansions can be 
worked out. For instance, one can analyze the structure of the differential cross section in 
order to discover in which observables D waves may show up. If one takes into account only 



S and P waves, then only To, 7\, and T2 coefficients contribute to Eq. (|2j). If we add D 
waves, two more quantities appear, T3 and T4. So the first place to look for D waves consists 
of checking to see if there is room for the appearance of these new terms. The currently 
available experimental data [9| can be described quite well using only T , T 1; and T 2 , and 
no T 3 or T 4 contribution appears to be required at present in the near-threshold region. 
Hence, any contribution of D waves to <Jt{s, 0) should appear only as a modification of To, 
T or T2. On the other hand, To and T2 are dominated by diagonal terms involving the 
multipoles (which can be immediately read from Tables [TTT1 and PVT) . namely |M 1+ | 2 , |M!_| 2 , 
and |i?o+| 2 , and thus any interference with D waves would be negligible compared with the 
leading-order terms. On the other hand, T% is entirely due to multipole interferences (all the 
coefficients in Table HVl are off-diagonal), and so any D-wscve interference with the dominant 
Mi+ multipole is a candidate for a non-negligible D-wave contribution to the observable. 
This affects any multipolar extraction from data using only S and P partial waves as was 
shown in Ref. 131 ]. 

It is very important to realize that without a so-called complete experiment — which 
at this point is not feasible in the near-threshold region — the multipole extraction from 
experimental data depends on the approach employed. For this reason we explore three pre- 
scriptions for the multipoles, with different phenomeno logical and theoretical content. The 
first one is a pure phenomenological approach where the energy dependence is prescribed 
for the multipoles (embedded in Sec. IIII Al) . The second consists in applying HBCHPT, 
which is the theoretically soundest approach to computing the S- and P-wave multipoles 
(Sec. IIIIBI) . The last one uses HBCHPT to compute the P waves, but a unitary prescrip- 
tion is employed for E 0+ (Sec. IIII Cj) — the reasons for this choice will be detailed in the 
corresponding section. 

For all three approaches the D waves are computed using the customary Born terms 
equivalent to the Born contribution to HBCHPT) and vector-meson exchange (u and p) 



181 ] . For the u and p parameters we have used those given by the dispersion analysis of the 



form factors in Ref. 19] that agrees with the latest analysis in Ref. [201 ] . The vector-meson 
correction is very small and the inclusion of D waves in this way is almost equivalent to the 
zeroth order in HBCHPT. For all approaches we perform fits using either solely S and P 
waves (SP fits) or S, P and D waves (SPD fits). The inclusion of the vector mesons in this 
fashion poses a model uncertainty in the calculation; however, this uncertainty is very small 
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and the conclusions on the effect of D waves in the observables and the extraction of the 
multipoles is effectively model independent. 

Within this framework we have performed fits to the latest experimental data from MAMI 
{^j] (171 differential cross sections and 7 photon asymmetries, spanning the energy range from 
threshold up to 166 MeV). As a fitting procedure we have used a hybrid optimization code 
:>ased on a genetic algorithm (GA) combined with the E04FCF routine from the NAG library 



2l| . In recent years, increased credit is being given in nuclear and particle physics to modern 



optimization procedures 22|, |23( and the error analysis for the parameters resulting from the 



fits. Modern and sophisticated optimization techniques such as GAs |24j have been developed 
over the past twenty years and have been applied to problems that were impossible to face 
with conventional tools. Although GAs are computationally more expensive than other 
algorithms, in a minimization problem they are much less likely to get stuck in local minima 
than are other approaches such as gradient-based minimization methods, and they allow one 
to explore a large parameter space more efficiently. Thus, in a multiparameter minimization 
as the one we face here, GAs probably provide the best possibility in searching for the 
minimum. Moreover, they provide additional information on the local minima structure. 
The details of the fitting procedure, its technicalities and advantages can be found in Ref. 



23]. 



III. THE ELECTROMAGNETIC MULTIPOLES 



A. Unitarity and the General Form of the Multipoles 



From general principles such as time reversal invariance and unitarity the 5" wave can be 
written as the combination of a smooth part and a cusp part 

E 0+ =e l5 ° [A + iPq+/m n+ ] ; W > W thr (n + n) 
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26] 



E 0+ =e l5 ° [A -P\q+\ /m K+ ] ; W < WW(vr + n) 



(6) 



where Sq is the 7i°p phase shift (which is very small), W is the invariant mass, W t hr{' 1T+n ) 
the invariant mass at the n + n threshold, q + is the 7r + center-of-mass momentum, A is 
E 0+ in the absence of the charge exchange re-scattering (smooth part), and j3 parameterizes 
the magnitude of the unitary cusp and can be calculated 25) on the basis of unitarity. 



Equations (jfJJ) provide a generalization of the Fermi-Watson theorem 



271 ] by removing the 



requirement of isospin conservation; they have been derived using a three-coupled-channel S- 



matrix approach in which unitarity and time reversal invariance are satisfied [25|, [26| . These 
calculations take the static isospin breaking (mass differences) as well as 7rN scattering to all 
orders into account. In the electromagnetic sector they have been carried out to first order 
in the fine structure constant a. The resulting equation for f3 is 

P = Eo+^p -> iT + n) x a(ir + n -> ir p) , (7) 

where a(7i + n — > 7i°p) is the ttN S'-wave charge exchange scattering length. The value of (5 
will be discussed in Sec. IIV Al 

The smooth part of Eq + can be approximated by using a Taylor expansion, and hence 
the entire S'-wave multipole can be written: 

p _ p(0) , F (l) k l ~ k l , -a Q*+ /ox 

Eo+ - E 0+ + E 0+ — - + — , (8) 

where E^J and E^ are free parameters, is the photon energy in laboratory frame, and 
= 144.681 MeV is the photon energy at threshold in the laboratory frame. The ir + 
center-of-mass momentum, is real above and imaginary below the 7r + threshold; this is 
a unitary cusp. HBCHPT matches this expansion for the S wave. 
Another expression for the E 0+ multipole is [jj, Q, [28] 



E 0+ = a + bJl-E*/E*, (9) 



where a and b are fitted to data, E n is the center-of-mass photon energy and E c the center- 
of-mass photon energy at the tt + production threshold. Equation is nearly equal to Eq. 
(JSJ) at the 3.5% level or better if b = (3 and E^_ = 0. Due to this latter restriction we do 
not use Eq. (Q in this work. 

The P waves can also be studied expanding in a Taylor series and assuming the multipoles 
are real we only need the lowest two orders: 

P i n(0) . D (i) k i ~ k i 



P U, + P 1)J1 ^; 2 = 1,2,3, (10) 

q 7T /m 7T o m n o 

where q w is the outgoing pion momentum and Pf^ and P^ are the coefficients of the 
expansion (free parameters). This expansion matches the energy dependence of HBCHPT 
for the real part of the multipoles. The Pi waves are related to the standard electromagnetic 
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multipoles through 



E 1+ = (P 1 + P 2 ) /6 (11) 
M 1+ = (P 1 - P 2 ) /6 + P 3 /3 (12) 
Mx_ = (P 3 + P 2 -Pi)/3 (13) 

Another choice for the P waves is the one made in Schmidt et al. 9|, where the prescrip- 
tion Pi = Pj e k^q^/m^o was employed where Pf is a constant fitted to data. We do not 
use such aprescription because it does not match the energy dependence of HBCHPT near 
threshold [3j and it is not truly the lowest-order Taylor expansion. The P-wave expansion is 
divided by q n due to the angular momentum barrier. The P-wave expansion, to first order, 
is 

A = D^ql/mlo , (14) 

where Pj = E 2+ , P2-, M 2+ , M 2 _, and q\ accounts for the correct angular momentum barrier. 
In this study the D waves are computed using Born terms and vector mesons that closely 
approximates this behaviour. 

Even fixing the D waves, these simple phenomenological prescriptions for the partial 
waves have too many parameters describing the energy dependence to perform a unique fit 
of the data. The reason is that the £ (polarized photon asymmetry) measurement has only 
been published at one energy. However, a more extensive asymmetry data set has been 
obtained in experiments recently performed at MAMI and the data analysis is in progress 
29I . Therefore, at this time a phenomenological analysis is not possible and one has to rely 
on theoretical approaches such as the ones we employ in the rest of the article. 

B. Heavy Baryon Chiral Perturbation Theory 

At present, the best available theoretical framework for studying pion photoproduction in 
the near-threshold region is HBCHPT. Because this approach is the one that was employed 



in 13j and the results presented here are the same as in that reference, we restrict ourselves 



to summarizing the approach and refer the reader to Ref. [13( for a discussion of the impact 
of D waves in the extraction of the E 0+ multipole, of the assessment of the LECs, as well as 
of the stability of P waves against the inclusion of D waves. In this section we summarize the 



10 



HBCHPT approach and leave the discussion of the impact of D waves in physical observables 
to Sec. IVl 

The explicit formulae for the S and P multipoles to one loop and up to 0(q*) can be found 
in Refs. |3|, |4| and constitute the starting point of our analysis. Due to the order-by-order 
renormalization process six LECs appear, and five have been fitted to pion photoproduction 
data: a\ and a 2 associated with the Eq + counterterm b p associated with the P3 multipole 
together with £ x and £ 2 associated with P-\ and P 2 , respectively. The c 4 LEC associated 
with Pi, P 2 , and P3 has been taken from [30] where it was determined from pion-nucleon 
scattering inside the Mandelstam triangle. Some other parameters appear in the calculation, 
but these are fixed. The full list is: the pion-nucleon coupling constant g n ^ — 13.1; the weak 
pion decay constant f n = 92.42 MeV, together with the anomalous magnetic moments of 
the proton and neutron, the nucleon axial charge g& (which we fix using the Goldberger- 
Trieman relation = g^Nf^/M^; and the masses of the particles. D waves are included 
as discussed in the previous section. 

In summary, using this approach, in Ref. 13j it was found that, contrary to what is 
customarily claimed in the literature, the SP approximation is not sufficient to obtain a 
complete description of the differential cross section in the near-threshold region and to 
extract the electromagnetic multipoles reliably. The inclusion of D waves does not affect 
the extraction of the P-wave multipoles but makes a significant difference where the E Q+ 
extraction is concerned, especially at and above the unitary cusp. The absence of D waves 
in the analysis affects both the Po+ multipole extraction and the determination of the 
associated LECs. 



C. Unitary Fit 

The Unitary fit is a hybrid approach that employs HBCHPT to compute the P waves 
and E 0+ using Eq. (jSJ) fitting the parameters Eq+ and Pq+ to the experimental data. These 
parameters have straightforward relationships to the LECs a% and a 2 or, more precisely, to 
their combinations a + = ai + a 2 and a_ == ay — a 2 . E^v fixes the value of the multipole 



at threshold and so does a + in HBCHPT [l3]. In the same way, Pq+ i s connected to the 
combination a + — a_ = 2a 2 , and thus one has a one-to-one relationship between both sets 
of parameters. 
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The HBCHPT approach has several shortcomings regarding the Eq + multipole. The first 
is its slow convergence and the second the lack of unitarity of the amplitude. In HBCHPT 
up to one loop and O (q*) the value of (3 = 2.71 x 10~ 3 /m n + is fixed by the imaginary part 
of E 0+ — which is parameter-free — and is not close to the unitary value of (3 that can be 
obtained from ttN scattering and whose value will be discussed in Sec. IIV Al The unitarity 
violation in HBCHPT is due to the truncation at the one loop (single rescattering) level. It 
is well known that Eq + does not converge very well ; however, as will be shown in Sec. 
IIVt the values for R,eE Q+ are very close for the HBCHPT and Unitary fits. 

The number of parameters is exactly the same as for the HBCHPT approach — five, 
namely, E^_ and E^J for the S wave and b p , £i and £2 for the P waves, which are the same 
as those defined in Sec. IIIIBI The D waves are included in the same way as in Sec. IIIIBI 
Under this approach we have performed two fits, one with D waves and another without 
them. The results are presented and discussed in Sec. HV] 



IV. RESULTS 



A. Multipole Extraction from the Experimental Data. 

In Table HT1 we summarize some results for the Unitary and HBCHPT approaches. All of 
the x 2 /DOF are about the same and are compatible within a 70% confidence level. In Fig. 
H the high quality of the SP Q and SPD fits for HBCHPT can be seen. The Unitary 
fits yield almost identical results. The P waves are quite alike from one fit to another with 
differences that lie below the 1.5% level at threshold. Moreover, the greatest differences are 
actually at threshold and when the energy increases the differences lie systematically below 
the 1% level, confirming the stability of P waves against the inclusion of D waves in the two 
approaches. In Table [IT] we also provide the values at threshold and at the unitary cusp for 
the real part of E 0+ . In the next paragraph and with the aid of Figs. [2] and [3] we discuss 
these values. 

The main difference between the Unitary and the HBCHPT approaches stems from the 
associated value of the unitary parameter (3 — see Eq. (J^J) — and how this affects the S'-wave 
extraction. Note that the sign of (3 is observable, not just its magnitude, and agrees with 
what is expected. The best experimental value of a(ir~p — > n°n) = —(0.122 ± 0.002)/m 7r +, 
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1 ^ — 1 1 1 1 1 

140 145 150 155 160 165 170 

E y (MeV) 

FIG. 1: Differential cross section for a fixed energy of £L = 159.9 MeV (upper figure) and a fixed 
angle of 6 = 90° (lower figure). Curve conventions: Solid: HBCHPT, SPD fit; Dashed: HBCHPT, 
SP fit. Both curves almost completely overlap. 

obtained from the observed width in the Is state of pionic hydrogen {3]]], was used. This is 
in excellent agreement with HBCHPT predictions of —(0.130 ± 0.006)/m 7r + }32| . Assuming 
isospin is a good symmetry, a(ir + n — > 7r°p) = —a(ir~p — > 7r°n). The latest measurement 
for Eo+frp -> 7r + n) = (28.06 ± 0.27 ± 0.45) x 10- 3 /m n + [33] (where the first uncertainty is 
statistical and the second is systematic) is in good agreement with the HBCHPT prediction 
of (28.2 ± 0.6) x 10 _3 /m 7r + [34]. These experimental values and the relationship given above 
lead to a value of /? = (3.43 ± 0.08) x lO -3 /™^ - 1 

1 Actually, it is energy-dependent, although for the current experimental accuracy it is a very good approx- 
imation to assume it is a constant. 
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TABLE II: Summary results for the Unitary and HBCHPT approaches. Units of Eq + in 10 3 /m^ 
and Pijq^ in lO^GeY" 1 /™-,^. 





HBCHPT fSP") 


HBCHPT fSPD"! 


Unitarv CSP > ) 


Unitarv (SPD) 


[ReKn , l 7r ° thr 


— 1 099 


— 1 0Q0 




— 1 196 


fRe/?n_ L r +thr 


—0 478 


—0 393 


—0 420 


—0 336 




73.88 


73.99 


72.99 


73.03 




-71.84 


-72.42 


-71.07 


-71.48 




76.71 


76.77 


75.98 


76.36 


X 2 /DOF (minimum) 


1.23 


1.25 


1.24 


1.23 


X 2 /DOF(70% C.L.) 


1.27 


1.28 


1.27 


1.26 


X 2 /DOF (90% C.L.) 


1.29 


1.30 


1.29 


1.28 



The results obtained in Ref. 13] using HBCHPT regarding the impact of D waves in 
the extraction of the Eq + multipole are confirmed using the Unitary approach in Sec. IIII CI 
D waves affect the extraction of the E 0+ multipole and if they are not included a reliable 
extraction cannot be achieved. When comparing the i?o+ extraction using HBCHPT and 
the Unitary approach it is better to split the analysis above and below the 7r + threshold. 
Regarding Rei?o+, from Fig. [2] it is evident that above and at the tt + threshold the extractions 
are very similar if the same number of partial waves is used, with only small differences that 
lie within the uncertainties. However, if one focuses on the region between the ir° and 7r + 
thresholds the situation is the opposite. As expected, due to the g 2 dependence — Eq. f fT4"|) 
— D waves have a negligible impact at 7r° threshold, as was shown in Ref. 13] for HBCHPT. 
As shown in Fig. [3l there is a sizeable difference between the two approaches regarding the 
imaginary part of Eq + , with the Unitary approach prediction being larger than that for 
the HBCHPT approach due to the larger value of (3. The piece of E 0+ that contributes 
to the imaginary part above the unitary cusp contributes to the real part below it. This 
means that this contribution in the real part of E 0+ near threshold is larger for the Unitary 
approach than for HBCHPT (this is the effect that we find at threshold) and makes the 
value of [Re.E'o+] 7r * hr 111 Table HT1 larger in absolute value for the Unitary approach than for 
HBCHPT. Hence, an accurate determination of the actual value of Eq + at threshold would 
help one to constrain the imaginary part of the multipole and vice versa. 
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FIG. 2: Extracted Rei?o+ multipole. Curve conventions: Thick Solid: HBCHPT, SPD fit; Thin 
solid: HBCHPT, SP fit; Thick dashed: Unitary, SPD fit; Thin dashed: Unitary, SP fit. 




0.14 0.145 0.15 0.155 0.16 0.165 0.17 

E y (GeV) 

FIG. 3: lmE 0+ multipole. Curve conventions: Solid: HBCHPT (SP and SPD fits); Dashed: 
Unitary (SP and SPD fits). 
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FIG. 4: Extracted Rei?o+ multipole ratio between the SP and SPD fits. Curve conventions: Solid: 
HBCHPT; Dashed: Unitary. 

In Fig. H] we focus on the effect of D waves, comparing the SP and the SPD fits for 
both HBCHPT and Unitary approach by computing the ratio for the real part of E 0+ . The 
impact of D waves is astonishing, changing the value of the multipole by almost 20% at the 
unitary cusp and reaching a 35% change at 165 MeV. 



In Ref . 



28] the Eq+ multipole [employing the phenomenological formula in Eq. was 



obtained fitting the subtraction constants to the dispersion analysis of Ref. [35J to compute 
the Fubini-Furlan-Rossetti (FFR) sum rule 36j] . This approach has the shortcoming that it 
relies on the extraction made through a dispersion analysis that is model dependent because 
it relies on MAID03 [37| to obtain the A\ amplitude that enters in the FFR sum rule. It does 
have the advantages that HBCHPT has better convergence inside the Mandelstam triangle 
and that all partial waves are incorporated in the dispersion analysis. 



B. Impact of D waves in Observables 

If the physical observables (differential cross section and asymmetries) are computed using 
either HBCHPT or the Unitary approaches, we do not find meaningful differences except 
for the case of TRO observables (T, P, G, and H; see, Figs. [5] and [6]) and, among them, 
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FIG. 5: T asymmetry for a fixed energy of E~ = 159.5 MeV (upper figure) and a fixed angle of 
9 = 90° (lower figure). Curve conventions: Solid: HBCHPT, SPD fit; Dashed: Unitary, SPD fit. 

the only one that is measurable in the near-threshold region is the T asymmetry 2 defined as 
the difference between the cross sections for a target polarized up and down along the axis 
perpendicular to the scattering plane. In Fig. Owe compare the results obtained for T using 
the Unitary and the HBCHPT approaches including D waves. These differences make the 
T asymmetry an excellent observable to test unitarity and to extract the imaginary part of 



E 0+ , because the latter is largely responsib 
to measure this asymmetry both at HI7S 



e of the difference. At present there are proposals 



loj and MAMI 



Henceforth we will focus the discussion on the effects of D waves on observables using 



2 Do not confuse the T asymmetry symbol with the Ti which stand for the partial wave expansion of the 
differential cross section. The latter has a subscript. 
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FIG. 6: TRO observables P, G, and H at outgoing pion angle 9 cm = 90° for SP and SPD fits using 
HBCHPT and the Unitary approaches. Curve conventions: Solid: SPD fit; Dashed: SP fit. Data 
are from [9J. In Fig. (b) both curves almost completely overlap. 



the HBCHPT calculation. This particular choice does not affect either the discussion or the 
conclusions. 

In Fig. [H]we present the comparison between the HBCHPT and Unitary approaches, both 
SP and SPD fits, for the TRO observables: recoil nucleon polarization asymmetry P, and 
double beam target polarized asymmetries G and H all as functions of the photon energy 
and for a fixed angle of 90° for HBCHPT and the Unit ary approach. The results are similar 
at other angles. The definition can be found in Ref. 15] or in Ref. [loj]; in this article we 
follow the sign conventions of Ref. [la ]. 

The effect of D waves in the G asymmetry is striking, but unmeasurable at present, 
since it is so small. The same can be said about P and H: the necessary accuracy to 
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measure the effect is at present beyond the current state-of-the-art. Differences between the 
Unitary and HBCHPT approaches are also unmeasurable at present. The most promising 
observables are the E and F asymmetries and in Fig. [7J we present them in terms of their 
energy dependence, which highlights the .D-wave effect for three different pion production 
angles: 45°, 90°, and 135°. When the SP and SPD fits yield different results, the predictions 
obtained using the HBCHPT and Unitary approaches are found to be significantly closer 
together than these differences. Therefore, in the remainder of the article we restrict the 
discussion to the HBCHPT fits. 

The combined measurement of the E and F asymmetries provides a good candidate for 
determining the D waves. The reason is that both are sizeable and behave in opposite ways 
for the D-wave inclusion: for those angles where D waves do not affect the E asymmetry 
they do impact the F asymmetry and vice versa (see Fig. [7]). More explicitly, at 90° the 
F asymmetry yields sizeable differences for the SP and SPD predictions, while not at 45° 
and 135°. The E asymmetry happens to be just the opposite. Hence, an analysis of both 
asymmetries using the SP and SPD approximations should be better at showing up the 
effects of D waves. 

In Fig. [8] we display the angular dependence of E and F for an intermediate energy 
of E 1 = 159.5 MeV. 3 The differences between results for the SP and SPD fits are not 
very apparent, but become more so when the observable is integrated over some convenient 
angular range (9i,9f); for instance, for the E asymmetry one might use 

f! f w E (8,e) sm ode 

<E>= J -^ — , (15) 

J°fW T (s,9) sin 9d9 

where Wt and We are defined by Eqs. (1A5I) and (1A7I) . This allows one to take advantage of 
cancellations in the asymmetry, thereby enhancing the features in which we are interested, 
and allows one to optimize the value of the figure of merit (FOM) — for the E asymmetry 
the FOM = ot < E > 2 — and hence minimize the statistical uncertainties. The observable 
< F > can be defined in a similar way for the F asymmetry by replacing We in Eq. (|15j) 
with Wf from Eq. (IA8I) . In Figs. [9] and [10] we provide < E > and < F > integrated over 
the appropriate angular ranges for the Crystal Ball at MAMI that maximize their FOMs 



3 We have chosen this energy both because it is intermediate and because it is the only one for which there 
are available £ data. 



19 



E F 




140 150 160 170 140 150 160 170 

E y (MeV) E y (MeV) 

FIG. 7: E and F asymmetries for SP and SPD fits using HBCHPT. Curve conventions: Solid: 
HBCHPT, SPD fit; Dashed: HBCHPT, SP fit. 




The main contribution to the differences between the SP and SPD fits in < E > are due 
first, to the differences in the E 0+ multipole and, second, to the interference of the M 1+ 
multipole with E2-. Again, the interference of D waves with M\ + constitutes the leading 
.D-wave correction. However, the fact that the biggest deviation is due to the modification 
of the E 0+ multipole proves that D waves are effectively entangled with the 5* wave and that 
the extraction of any of them must rely on the knowledge of the others. 

These asymmetries are large and actually would allow one to find a trace of D waves. 
However, we believe that the combined measurement of both would improve the analysis 
beyond doubt. Nevertheless, the angular distribution is also of great interest in order to 
obtain the coefficients Fi and Ei in Eqs. ( IA8I) and (1A7j) . There are plans to measure F at 
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FIG. 8: Angular dependence for the E and F asymmetries at E~ = 159.5 MeV. Curve conventions: 
Solid: HBCHPT, SPD fit; Dashed: HBCHPT, SP fit. 



HI 7 S [10J and MAMI [38 1. 

At first thought, the photon beam asymmetry £ seems to be a promising observable to 
pin down D waves if enough accuracy is achieved. If we consider only S and P waves, this 
asymmetry can be written: 

- WtE = S sin 2 9 , (16) 

where So can be read from Table IVIIIl The inclusion of D waves modifies So and adds 
two new terms Si V\ {&) sin 2 9 and S2V2 (0) sin 2 8 to the right-hand side of Eq. (fIBl) . The 
leading term would be Si; however, this term seems to be very small and difficult to isolate, 
as we show in Fig. [TT] where we compare the SP and SPD fits with experimental data 
Moreover, the D-w&ve contribution could be easily hidden by the fitting, embedding the D 
waves in the effective parameters of the multipoles. In order to avoid this situation, some 
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FIG. 9: E asymmetry depending on energy and integrated in the angular range from 90° to 140°. 
Curve conventions: Solid HBCHPT, SPD fit; Dashed: HBCHPT, SP fit. Errors are computed at 
the 70% confidence level. 
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FIG. 10: F asymmetry depending on energy and integrated in the angular range from 70° to 140°. 
Curve conventions: Solid: HBCHPT, SPD fit; Dashed: HBCHPT, SP fit. Errors are computed at 
the 70% confidence level. 
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FIG. 11: X asymmetry for a fixed energy of E~ = 159.5 MeV (upper figure) and a fixed angle of 
6 = 90° (lower figure). Curve conventions: Solid HBCHPT, SPD fit; Dashed: HBCHPT, SP fit. 

other observables should be measured in order to find inconsistencies in the SP fits that can 
only be solved by adding new terms (D waves). 

For this reason, one should also consider some other convenient observables. For instance, 
if the E asymmetry is measured together with the differential cross section and enough 
angular points are acquired, it is possible to extract the values of T, and Ei at a certain 
energy. If we compare the tables for T\ and E% (Tables IIVI and Kill) we notice that both 
have the same coefficients for the S-P interference, and hence any difference between them 
is due to higher partial waves (D waves). Moreover, if we assume E\ + ~ (which is actually 
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FIG. 12: OS = T 1 -E 1 for HBCHPT SPD fit (solid) and HBCHPT SP (dashed). The SP calculation 
yields zero and is shown for comparison purpouse. 

a sensible assumption), it is straightforward to obtain 

23 =Ti - E l 

18 ( 17 ) 
~— Re{M* + [4 (M 2+ - £ 2+ ) + M 2 _ + £ 2 _]} , 

whose measurement provides a perfect test of the contribution of D waves because the P 
wave M\ + is large and well known. 

Both pieces of information can be measured and one can obtain an unmistakable trace of 
D waves in low-energy pion photoproduction due to the fact that under the SP approxima- 
tion this difference is exactly zero. Hence, any deviation from zero will be due to the interfer- 
ence of P and D waves. Moreover, as long as M\ + is the dominant P-wave contribution (and 
well known) and E i+ is negligible, the D-w&ve combination 4 (M 2+ — E 2 +) +M 2 _ +E 2 - can 
be isolated, providing a test of its presence or absence. In Fig. [12] we show our prediction for 
23, whose magnitude is possibly measurable with the current state-of-the-art facilities. This 
would also provide insight on the D waves, testing how accurate is the Born+vector-meson 
approach employed in this article. 

It is worth noticing an important difference between this observable and the measurement 
of E and F in the search of D waves. This combination is truly a measurement of D waves, 
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FIG. 13: Ex for fits HBCHPT SPD (solid), HBCHPT SP (dashed), and HBCHPT SPD (dotted) 
with the D waves removed. 

while in E and F what is involved is the impact of D waves on the extraction of the S- 
wave. Actually, the main contribution to the difference between the SP and SPD fits is 
not the addition of D waves, it is the modification of the S wave due to the addition of D 
waves and the subsequent impact of this difference on the observables. Figure [13] illustrates 
this statement. We display the results from three calculations for the E\ observable. The 
solid line is the HBCHPT SPD fit and the dotted one the HBCHPT SP one. The difference 
between the two is apparent. The third curve is the HBCHPT SPD fit without the D waves. 
As long as P waves are almost equal for the two calculations, the difference between the 
dotted and dashed curves is due to the difference between the S waves, and the difference 
between the solid and the dotted curves is the true contribution of the D waves, mainly 
through their interference with P waves. 

Another good candidate for extracting the E 0+ multipole is 

£ = T + 2E = 3\E 0+ \ 2 + 3|Mx_| 2 

(18) 

+ 12|^ 1+ | 2 + 12Re{£* + M 1+ } + 5D , 

where 5D stands for the .D-wave content, namely D waves interfering with D waves, which 
guarantees a very small contribution to the observable. As long as the P waves are well 
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known from HBCHPT and E\ + is very small, then we can extract the value of |Po+| 2 if we 
know the P-wave Mi_. If we approximate D ~ 3 1 1 2 + 3|Mi_| 2 , the error is less than a 
6%. Combined with other measurements, D would help to pin down the E 0+ multipole. 

Actually, the particular choice 2) can be exploited in a different way Experimentally, 
asymmetries can be measured more accurately than cross sections, and the quantities E , 
Ei, Fo and F\ can be measured with relatively small systematic errors. Moreover, we can 
define 



with Ti = Ti/Tq and T 2 = T2/T0 that can both also be measured very accurately. Hence, 
the key quantity to improve is T , which is the total cross section up to a kinematical factor 



From accurate measurements of the asymmetries (T, F and E) it is possible to extract 
the multipoles and use this knowledge to isolate T in Eq. (Tl8l) . providing a procedure 
to determine the cross section more accurately. These are a few of the analyses that a 
deep knowledge of the structure of the observables in terms of the multipoles allows one to 
undertake using the tables in Appendix [A] 

C. Contribution of Higher Partial Waves to the Observables 

At this point a fair question to ask is why should we stop at D waves and not include 
F waves in the analysis. Within the current experimental accuracy (and probably for years 
to come) the answer is that F waves do not play a role and can be dismissed. The reason 
to make this statement relies on the structure of the observables presented in Appendix [A] 
Looking at the expansions in Eqs. (IA5I) . (IA6I) . flA7j) and (IA8I) . coefficients such as T 3 , T 4 , 
S 2 , E 3 , F 4 , F 2 , and F 3 are negligible experimentally. Coefficients such as T , T 2 , So, E Q , 
and F\ are |Mi + | 2 and |Mi_| 2 dominated, and so if F waves would have impact the only 
observables that remain are those where D waves already make an impact, namely T\, S\, 
Ei, and Fo- At this point we remind the reader that the enhancement of D waves is due to 
their interference with the large P-wave multipoles, M\ + (primarily) and Mi_, and therefore 
any non-negligible F wave contribution should be related to the same kind of interference. 
Due to symmetry reasons, in T 1; E\, and F F waves interfere only with D and G waves. 



W T = T 1 + T{P X (9) + T 2 P 2 (#) + ... 



(19) 




(20) 
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The argument is as follows: first, for Ti, E\ and Fq only partial waves with opposite parity 
can interfere (see Tables HVt Kill and IXVIj) ; second, these partial waves can differ only by 
one unit of angular momentum. Si is the only coefficient where F waves interfere with P 
waves, but we have found that P-wave/P-wave interferences for this observable are small 
and hence we expect F-wave/P-wave interferences to be even smaller. So F waves are highly 
suppressed in the observables and can be safely neglected in the near-threshold region within 
the current experimental state of the art. Clearly our expectation is that G waves and higher 
play an even more minor role in the near-threshold region. 



V. SUMMARY AND CONCLUSIONS 

In Ref. 



13| neutral pion photoproduction from the proton was studied in the region from 
threshold up to 167 MeV using HBCHPT. The E 0+ multipole and D waves are entangled 
in such a way that the accurate determination of the former (including the extracted values 
of the LECs a\ and t^) relies on the proper inclusion of D waves. In this article we comple- 
ment the analysis in Ref. [l^l by employing a phenomenological approach that needs more 
experimental input to draw reliable conclusions and a Unitary approach that uses HBCHPT 
for P waves and a unitary prescription for the E 0+ multipole. 

We have found no impact of D waves at 7r° threshold on the S wave, although a noticeable 
difference is found between the Unitary and HBCHPT approaches due to the difference in 
the value of (3. We have also considered partial waves higher than L = 2 and have found 
that they are not important and can be neglected. We have found that P waves are stable 
at the 1% level against the inclusion of D waves or the use of a different prescription for the 
E 0+ multipole. 

The major difference between the SP and SPD calculations of the observables does not 
rely on the direct contribution of D waves, but on the modification of the Fo+ multipole, 
which is misleadingly extracted from data if only S and P waves are included. Although 
D waves are typically weak, the modification of the S wave is sufficiently large to alter 
the prediction for the double-polarization observables E and F enough to to be able to 
distinguish between approaches with and without D waves, leading to measurable effects 
that will allow one to pin down D waves experimentally. In the next generation of photo- 

these could be measured at MAMI and HI7S. Specifically, it 



pion experiments [10, 12 
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has been demonstrated that the statistical accuracy of the recently proposed polarization 
experiment to determine F at MAMI will be sufficient to unveil the importance of the D- 
wave contribution. In this article we have gone further and shown also that the E observable 
is sensitive to the .D-wave contribution and that the combination Ti — Ei will provide a direct 
quantitative measure of the D-wave effect (Fig. [T2"|) . We thus expect in the next few years 
a direct experimental test of the ideas put forward in this article. 

We also explored the influence of D waves in 7rN scattering and found no sizeable impact 
on the partial-wave extraction. In addition, the treatment of the Coulomb interaction poses 
uncertainties that are larger than the estimated .D-wave effects. 

In conclusion, D waves cannot be dismissed in the analysis of low-energy neutral pion 
photoproduction from the proton due to the soft nature of the S wave that is direct conse- 
quence of chiral symmetry and the Nambu-Goldstone nature of the pion. 
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APPENDIX A: PHYSICAL OBSERVABLES IN TERMS OF THE ELECTRO- 
MAGNETIC RESPONSES: MULTIPOLAR TABLES UP TO D WAVES 

The differential cross section and asymmetries can be written in terms of electromagnetic 
responses 



<tt(s,0) 



W t (s,8) 



(Al) 



E( s ,e) 



F(s,e) 



W s (s,0) 
W T (s,8) 
W E (s,0) 

w T (s,e) 

W F (s,0) 
W T (s,6) ' 



(A4) 



(A2) 



(A3) 
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The responses Wt, Ws, We, and Wf are defined in term of the electromagnetic multipoles 
up to D waves: 



W T = T (s) + 7\ (s) Vx (9) + T 2 (s) V 2 (9) 
+T 3 (s)P 3 (9) + T 4 (s)V 4 (9) 

W s (s,9) = [S (s)+S 1 (s)Vi (9) 

+ S 2 (s) V 2 (9) } sin 2 9 

W E = E (s) + Ex (s) V x (9) + E 2 (s) V 2 (9) 
+ E 3 (s)V 3 (9) + E 4 (s)T 4 (9) 

W F (s,9)=[F (s) + F 1 (s)V 1 (9) 

+ F 2 (s)V 2 (9) +F 3 (s)V 3 (9) ]sin0 
where Pj (9) are the Legendre polynomials in terms of cos 9 and 

T n (s) = "£Re{M*(s) T l J Mj (s) } 

ij 

S n (s) = J2MM*(s) SUM^s)} 

ij 

E n (s) = J2Re{M*(s) BftMiia)} 

ij 

F n (s) = "£Re{M*(s) F%> Mj (s) } 



ij 



(A5) 



(A6) 



(A7) 



(A8) 



(A9) 
(A10) 
(All) 
(A12) 



where Mj (s) = E 0+ , E 1+ , E 2+ , E 2 _, M 1+ , M x _, M 2+ , M 2 _. 

The coefficients T% are provided in Tables EH US El EH and EEQ S% in Tables IViTTI 
US and S E% in Tables Kfl ED IXTTT1 IXT71 and EH and F^' in Tables IXVTl IXVHl IXVilil 
and IXIXl The equations necessary to compute the tables can be found in Ref. |l4|. For a 
complete list of polarization observables we refer the reader to Refs. [ij], fl6 [- 
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